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ABSTRACT 

The BESS-Polar Collaboration measured the energy spectra of cosmic-ray protons and helium 
during two long-duration balloon flights over Antarctica in December 2004 and December 2007, at 
substantially different levels of solar modulation. Proton and helium spectra probe the origin and 
propagation history of cosmic rays in the galaxy, and are essential to calculations of the expected 
spectra of cosmic-ray antiprotons, positrons, and electrons from interactions of primary cosmic- 
ray nuclei with the interstellar gas, and to calculations of atmospheric muons and neutrinos. We 
report absolute spectra at the top of the atmosphere for cosmic-ray protons in the kinetic energy 
range 0.2-160 GeV and helium nuclei 0.15-80 GeV/nucleon. The corresponding magnetic rigidity 
ranges are 0.6-160 GV for protons and 1.1-160 GV for helium. These spectra are compared to 
measurements from previous BESS flights and from ATlC-2, PAMELA, and AMS-02. We also 
report the ratio of the proton and helium fluxes from 1.1 GV to 160 GV and compare to ratios 
from PAMELA and AMS-02. 
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1. Introduction 

We report new measurements of the energy 
spectra of cosmic-ray protons and helium made by 
the BESS-Polar Gollaboration (Balloon-borne Ex- 
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periment with a Superconducting Spectrometer - 
Polar) during long-duration balloon (LDB) flights 
over Antarctica in December 2004 (BESS-Polar I), 
prior to the last solar minimum, and in Decem¬ 
ber 2007 (BESS-Polar II), at solar minimum. The 
trajectories of both flights are shown in Figure [T] 
The absolute fluxes, corrected to the top of the at¬ 
mosphere (TOA), are reported in energy ranges of 
0.2-160 GeV for protons and 0.15-80 GeV/nucleon 
for helium. In magnetic rigidity (i? = pc/.Ze, 
momentum divided by electric charge), the cor¬ 
responding ranges are 0.6-160 GV for protons and 
1.1-160 GV for helium. These new spectra are 
compared to results previously repor ted by BESS 


and t o measurem ents from ATIC-2 (Panov et ah. 
2007ll . P AMELA jAdriani et al.]l201lL l20l4 , and 


AMS-02 dAguilar et al.||2015alM . We combine our 
measurements to give the ratio of the proton and 
helium fluxes from 1.1 GV to 160 GV and compare 
this to ratios reported by PAMELA and AMS-02. 


The BESS-Polar measurements reported here 
are very important to a wide range of astrophys- 
ical studies and calculations. Galactic cosmic 
rays (GGR) near the peak of the all-particle en¬ 
ergy spectrum at Earth, ^2 GeV, are made up 
of ^98% atomic nuclei (including antiprotons) to¬ 
gether with ^2% electrons and positrons, and the 
fraction of atomic nuclei increases rapidly with in¬ 
creasing energy. At energies where direct measure¬ 
ments are possible, to about the spectral break 
known as the “knee” at ^3 PeV, GGR nuclei 
are ^ 87 % protons and ~12% helium. The abso¬ 
lute fluxes and spectral shapes of primary cosmic- 
ray species are central to understanding the ori¬ 
gin and the propagation history of cosmic-rays in 
the galaxy and provide crucial tests of cosmic-ray 
models. Proton and helium spectra are essential 
inputs to calculations of the spectra of cosmic-ray 
antiprotons, positrons and electrons produced as 
secondary products of primary GGR interactions 
with the interstellar gas. They are also required 
for calculations of the fluxes of atmospheric muons 
and neutrinos that are important backgrounds 
to ground-based neutrino experiments and result 
from a decay chain beginning with the production 
of charged pions by interactions of GGR and at¬ 
mospheric nuclei. 



Fig. 1.— Flight trajectories over Antarctica of 
BESS-Polar I in 2004 (green) and BESS-Polar II 
in 2007/2008 (first orbit blue, second orbit red). 


2. The BESS-Polar Program 


The BESS ins trument ( Aiima et al.l I 2 OOOI : 
Haino et"^ l2004bll was developed as a high- 
resolution, high-geometric-acceptance magnetic- 
rigidity spectrometer for sensitive measurements 
of cosmic-ray antiparticles, searches for antinuclei, 
and precise measurements of the absolute fluxes 
of light GGR elements and isotopes. Between 
1993 and 2002, versions of the original instrument 
incorporating continuous improvements in perfor¬ 
mance flew 8 times from Lynn Lake, Manitoba, 
Ganada and once from Fort Sumner, New Mexico, 


USA (jMitchell et al.lDOOdll . 


BESS-Polar (Yamamoto et al.l2002l Yoshida et al. 

2004 

: Yamamoto & MitchelJ2013l: Yamamoto et al.l 

2013 

: Mitchell et al. 2013f) is the most advanced 


version, based on the general BESS instrument 
architecture but developed specifically for LDB 
flights over Antarctica to provide much higher 
statistics at consistently low geomagnetic cutoff 
than could be obtained in conventional balloon 
flights. BESS-Polar presents 45% less material to 
incident particles penetrating the full instrument 
because it has no outer pressure vessel, a thinner 
magnet coil, and thinner detectors. For the mea¬ 
surements reported here, this lowers the thresh¬ 
old energy to ^200 MeV for protons and ~150 
MeV/nucleon for helium. Completely revised low- 
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power electronics support extended LDB flights 
powered by photovoltaic (PV) arrays and measure 
the abundant protons and helium nuclei without 
the event selection required in BESS to control 
trigger rates. BESS-Polar I was launched from 
Williams Field near McMurdo Station, Antarc¬ 
tica, on December 13th, 2004 and flew for over 8.5 
days with the magnet energized at 37 km to 39.5 
km (average atmospheric overburden 4.33 g/cm^) 
and geomagnetic cutoff R below 0.2 GV, record¬ 
ing 2.1 4 terabytes of da ta on 9 x 10® cosmic-ray 


events ( Abe et al. 2008l l. 


BESS-Polar 11 incorporated considerable im¬ 
provements in instrument and payload systems, 
including improved magnet cryogen hold time, re¬ 
duced deadtime per event, and much greater on¬ 
board storage capacity. It was launched on De¬ 
cember 23, 2007 and observed cosmic rays for 24.5 
days with the magnet energized at 34 km to 38 
km (average of 5.81 g/cm^) and cutoff R below 
0.5 GV, accumu lating 13.6 terabytes of data on 
4.7 X 10® events ( Mitchell et al. 2013 1. 

Previous papers from the BESS-Polar program 
have reported precise m easurements of antiprotons 
( Abe et ahlbOOSl 201^ Jand a sensitive antihelium 
search ( Abe et ^ 2012b[l and discuss instrumen¬ 
tation specific to those measurements. 


3. The BESS-Polar Instruments 

In both BESS-Polar instruments (Figure [2|), a 
uniform magnetic field is produced by a thin su¬ 
perconducting solenoid magnet 0.9 m in diameter 
and 1.4 m long, incorporating a 0.8 m diameter, 
1.4 m long coaxial warm bore in which detectors 
are located. The magnet is capable of stable op¬ 
eration at 1 Tesla (1 T), but in both flights was 
operated at 0.8 T for additional stability margin. 
Drift-chamber tracking detectors in the warm bore 
measure R. The cylindrical configuration of the 
BESS-Polar spectrometer gives it a large angular 
acceptance. With ^1 m® time-of-flight (TOE) ho- 
doscopes separated by ~1.6 m, the nominal geo¬ 
metric acceptance of BESS-Polar is ^0.3 m^sr at 
10 GV. 

The primary TOE system composed of upper 
(UTOF) and lower (LTOF) plastic scintillator ho- 
doscopes measures velocity (/?) and differential en¬ 
ergy loss {dE/dx) and forms the event trigger. A 
thin plastic scintillator middle-TOF (MTOF) ho- 


doscope is installed on the lower surface of the 
magnet bore to measure low-energy particles that 
cannot reach the LTOF. The MTOF further low¬ 
ers the threshold energy to about 100 MeV for 
antiproton or proton measurements, but has not 
been used here because the low energy spectra are 
strongly affected by solar modulation and our fo¬ 
cus is on higher energies. 

A threshold-type silica-aerogel Cherenkov counter 
(ACC), below the spectrometer, helps reject e~ 
and backgrounds to the antiproton measure¬ 
ments. This is not used in the proton and helium 
analysis reported here because contamination by 
light backgrounds is insignificant. Its use in the 
antiproton analysis has been described previously 
dAbe et al.ll20nal2012ah . 

Particle electric charge (Ze) is determined using 
dE/dx and /3. Particle mass (m) is calculated from 
A, Z, and /3 as: 

m = i?Zev'l//32 - 1. (1) 


Charge sign is determined by the direction of cur¬ 
vature in the magnetic field. 


The superconducting magnet uses NbTi/Cu 
wire with high-strength aluminum stabilizer in- 
corporati ng Ni doped filame nts to increase its 
strength ( Makida et al.l[^09l l. The magnet and 
cryostat have only 2.46 g/cm^/wall. The mag¬ 
net operates in persistent mode with a cryogen 
hold time of 11 days in BESS-Polar 1, using a 400- 
liter liquid helium (LHe) reservoir, improved to 25 
days in BESS-Polar 11 with 520 liters of LHe. The 
tracking system uses a JET-type drift chamber 
(JET) with essentially no material in the track¬ 
ing volume, other than gas, to minimize scatter¬ 
ing. Inner drift chambers (IDCs) with two track¬ 
ing layers each are located above and below the 
JET. 


The JET and IDCs accurately determine the 
position of particles passing through each drift cell 
from the recorded drift times and measure their 
dE/dx from the digitized signal amplitudes at the 
sense wires. Pure CO 2 fills the drift chambers 
and fresh gas is circulated during flight by a semi¬ 
active flow control system and exhausted through 
the pressurized warm-bore of the magnet to en¬ 
sure that the drift gas is not contaminated. The 
signals of the JET and IDCs are read out using 
flash analog-to-digital converters (FADC) capable 
of recording both the time of detection at the sense 
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Fig. 2.— Crosssection and side views of the BESS-Polar II instrument 


wire and the shape and amplitude of the detected 
pulse. Particle tracks are measured by up to 52 
points (48 from the JET and 2 from each of the 
2 IDCs) with ^125 /xm resolution in the bend¬ 
ing {T(f)) direction (perpendicular to the magnetic 
field). The uniform magnetic field of the BESS- 
Polar solenoid makes the R measurements rela¬ 
tively insensitive to tracking accuracy in the 2 co¬ 
ordinate (parallel to the magnetic field), and this 
is accomplished using vernier pads in each layer 
of the IDCs, giving ^650 /xm resolution, and JET 
measurements using charge division to give ^24 
mm resolution. The continuous tracking with min¬ 
imal scattering provided by its tracking system 
enables BESS-Polar to readily distinguish inter¬ 
actions and hard scatters. A truncated mean of 
the integrated charges measured by the JET pro¬ 
vides an energy loss {dE/dx) measurements with 
a resolution of 10%. The maximum detectable 
rigidity (MDR, the highest R curved track distin¬ 
guishable from a straight line) is 240 GV for the 
longest tracks through the spectrometer, expect¬ 
ing 52 points to be used in fitting tracks. 

The scintillator hodoscopes consist of 10 UTOE 
and 12 LTOE Eljen E204 scintillator paddles, each 
0.95 m long with cross sections 96.5 mm wide 
X 10 mm thick in BESS-Polar I and 100 mm 
X 12 mm in BESS-Polar II. Hamamatsu R6504 
magnetic-field-tolerant, fine-mesh photomultiplier 
tubes (PMTs) at each end view the scintillators 
through UV-transmitting acrylic light guides. The 
timing resolution is 120 ps, giving a j3~^ resolution 
of 2.5%. The nominal BESS-Polar event trigger 
is a four-fold coincidence that requires a PMT on 
each end of both the UTOE and LTOF to exceed a 
fixed signal threshold. Coincidences of the UTOF 


and MTOF also trigger events. 

TOF signals are digitized by low-power common- 
stop time-to-digital converters (TDCs) attached 
to the PMT anodes and 16-bit charge-to-digital 
converters (QDCs) attached to both the 13th and 
18th dynodes to increase dynamic range. The 
18th dynode measurements are normally used 
to measure dE/dx for light particles and nuclei. 
As noted, JET and IDC signals are digitized by 
FADCs. All readout electronics are highly parallel 
and the full processing and communication time 
for each event is ^40 fis. Event data are recorded 
to onboard hard disk drives with a total capac¬ 
ity of 3.5 terabytes (TB) in BESS-Polar I and 16 
TB in BESS-Polar II. The data acquisition event 
rate was ^1.4 kHz in BESS-Polar I (due to PMT 
issues that reduced the acceptance, discussed be¬ 
low) and ~2.5 kHz in BESS Polar H. Total power 
consumption was ^420 W. Because the magnet 
fixed the orientation of the payload with respect 
to the local geomagnetic field, all sides of the pay- 
load were exposed to the sun each day and an 
omnidirectional photovoltaic array was required. 

4. In-Flight Instrument Issues 

The UTOF and LTOF operate at ambient 
pressure (~5 mbar) and protection against high- 
voltage (HV) breakdown was required. In BESS- 
Polar I this was accomplished by potting the high- 
voltage components. In the BESS-Polar I flight, 
18 of 44 PMTs in the TOE had to be turned 
off because their potting failed under a combina¬ 
tion of cold and vacuum, and they began to draw 
excessive current. This required the trigger to be 
changed in flight to a two-fold coincidence between 
any UTOF PMT and any LTOF PMT. Although 
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some paddles did not have a working PMT and 
did not provide trigger signals, this retained 74% 
of the nominal trigger geometric acceptance. In 
BESS-Polar II there were no HV issues due to 
vacuum, but 2 of the 44 PMTs had HV control 
problems and were turned off, one on the UTOF 
and one on the LTOF. With a change to a two¬ 
fold coincidence, all TOF paddles provided trigger 
signals and 100% of the nominal trigger geometric 
acceptance was retained. 

In BFSS-Polar I, the tracking system was stable 
throughout the flight. In BESS-Polar II, however, 
the JET HV suddenly became unstable after more 
than a day at float and the voltage and the cur¬ 
rent began to fluctuate. The HV was immediately 
reduced to 90% of the nominal operating value 
and the CO 2 pressure was reduced to keep the 
same drift velocity. This intervention succeeded 
and the chamber continued to function well for 
the remainder of the flight, although the HV fluc¬ 
tuated for periods between two distinct voltages. 
Despite this, normal tracking resolution was ob¬ 
tained for more than 90% of the science observa¬ 
tion time by development of algorithms that cal¬ 
ibrate the tracker over short time intervals and 
depend on its HV state. Time lost when the HV 
was unstable is reflected in the trigger live time 
discussed in Section [531 Straight tracks obtained 
in flight without magnetic field after the LHe was 
expended were used as an absolute reference of 
events with “infinite i?” to check the calibration 
and ensure that it was not biased. 


5. Data Analysis 


5.1. Tracking System Calibration 

Each track point in the drift chambers was cal¬ 
culated from the measured drift times in the JET 
drift cells and the IDCs. The position recon¬ 
struction started with space-time templates cal¬ 
culated by applying th e CERN GARFIELD simu¬ 
lation code ( Verrnholll998ll to the nominal design. 
These templates were adjusted to match the actual 
chambers by iteratively modifying the parameters. 
The calibration process considered FADC time 
offsets, electric-field and CO 2 pressure-dependent 
drift velocities, Lorentz-angle corrections, wire po¬ 
sition displacements, and field distortion correc¬ 
tions. Tracks through the magnetic field were con¬ 
structed from the measured points and used to 


calculate R. 

The limitations to the accuracy of the R mea¬ 
surements arise from the characteristics of the 
spectrometer including the integral magnetic field, 
field uniformity and absolute value, track length in 
the field, and tracking resolution. There are also 
potential uncertainties from deflection offsets due 
to small misalignments of the sense wires. 

Prior to flight, the magnet was energized to a 
central magnetic field of 0.8 T, and then put into 
persistent current mode. The subsequent decay of 
the magnetic field is very slow and it was essen¬ 
tially constant for the duration of the flight. Thus, 
the uncertainty in the absolute magnetic field is 
dominated by the uncertainty in the excitation 
current, estimated as 0.5%. The resulting system¬ 
atic uncertainties in the observed fluxes were esti¬ 
mated by varying the magnetic field ±0.5% from 
the nominal value and found to be < 1.5% for both 
protons and helium. 

The measured distribution of the deflection res¬ 
olution {a{R~^)) in both BESS-Polar I and BESS- 
Polar H has mean of 4.2 TV“^ for tracks with 
the longest possible lengths in the spectrometer. 
Over the full fiducial volume used in this analysis, 
the mean of the measured deflection distribution 
is 5.0 TV“^ corresponding to an MDR of 200 GV. 
The precise alignment of the JET with respect to 
the IDCs was determined by checking the consis¬ 
tency between tracks reconstructed by JET and 
points measured by IDCs. Wire position displace¬ 
ments were determined using ground muon data 
for BESS-Polar I and with data taken without 
magnetic-field at the end of the flight for BESS- 
Polar H. A small net displacement was observed 
in both cases. The wire position calibration of the 
JET gives a maximum sense-wire alignment uncer¬ 
tainty of <10 /xm. It was expected that this would 
be random and would have no significant effect on 
performance. However, a coherent misalignment 
of all wires by the maximum 10 /xm could intro¬ 
duce a systematic offset in the deflection measure¬ 
ment of =0.33 TV“^ and result in the ob¬ 

served deflection displacement. This was included 
as part of the systematic error attributed to cham¬ 
ber alignment. A detailed model of the struc¬ 
ture, operating conditions, and response function 
of the drift chamber tracker was implemented in a 
Geant4 Monte Carlo simulation and the simulated 
R resolution accurately reproduced the measured 
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resolution. 

5.2. Track Reconstruction and Accep¬ 
tance 

Tracks were fit by an iterative method devel¬ 
oped to give high track-finding efficiency with lit¬ 
tle rigidity dependence. An initial track was first 
generated using all valid chamber hits. This was 
then refined by removing from the htting all wires 
whose measured track position was more than 20 
artj, away from the corresponding value from the 
previous track fit and then re-fitting the track us¬ 
ing only the surviving points. This process was 
repeated for 15, 10, and finally 5 ar^ and the final 
track was fit using points from those wires that 
survived all cuts. 

For the best R resolution consistent with max¬ 
imum efficiency and geometric acceptance, parti¬ 
cles are required to traverse the tracker completely 
within a fiducial volume defined by the central six 
columns of drift cells (out of eight) in the JET and 
a requirement on the cosine of the zenith angle of 
cos^zenith > 0.80. This reduced the geometric ac¬ 
ceptance at 10 GV to 0.149 m^sr in BESS-Polar I 
(due to TOE PMT issues discussed in Section |4l) 
and 0.199 m^sr in BESS-Polar II. For both instru¬ 
ments this requirement gives >40 points expected 
from the JET and 4 points required from the IDCs. 
The error comes from the uncertainty of the de¬ 
tector alignment of less than 1 mm, corresponding 
to about 0.2% systematic error of the acceptance 
by UTOF and LTOF. 

5.3. Monte Carlo Simulated Efficiencies 

Selection criteria were applied as cuts to the 
event data. These are characterized by cut selec¬ 
tion efficiencies defined as the probability of an 
event surviving a particular cut. There are two 
classes of selection. Selections whose effects can¬ 
not be determined directly from the flight data 
are discussed here. Their corresponding selection 
efficiencies are estimated using Geant4 (Version 
lO.OO.pOl) Monte Carlo simulations and are com¬ 
bined into an overall selection efficiency ffMC- Se¬ 
lections whose effects can be directly determined 
from flight data are discussed in Section [5^ 

The selection criteria were almost identical in 
BESS-Polar I and BESS-Polar II analyses with the 
following exceptions. In the BESS-Polar I analy¬ 


sis, we allowed only one hit paddle on either the 
UTOF or LTOF to be read out by a single PMT 
and required the hit paddle on the other TOF layer 
to be read out by 2 PMTs. The corresponding se¬ 
lection efficiency was 59% for protons and 62% for 
helium. In BESS-Polar II, 2 of the 34 IDG wires 
were excluded because of low signals due to the 
HV supply issue. To ensure the best TOF reso¬ 
lution and to allow z position constancy between 
the TOF and tracking to be used, we conserva¬ 
tively required each TOF paddle to be read out by 
2 PMTs. The selection efficiency from the combi¬ 
nation of these cuts were 58% for proton and 61% 
for helium. 

To eliminate rare events in which an interaction 
took place in or near the instrument, single-track 
events were required having only one isolated track 
and no more than two paddles hit in each layer of 
the TOF hodoscopes. The non-interaction selec¬ 
tion efficiency as a function of R was estimated by 
applying the single-track criteria to events from 
Geant4 simulations. At 10 GeV/nucleon it was 
about 90% for protons and 77% for helium. The 
error in the non-interaction efficiency was esti¬ 
mated by comparing the number of LTOF paddles 
hit in the flight data and in the simulated data. 

The resolution of the IDGs deteriorates for 
tracks passing both near the sense wires and at 
large drift distances. To ensure the best tracking, 
an effective IDG fiducial region of about ^85% of 
the active area was defined by cuts on the drift 
time to exclude the regions very near the sense 
wires (~10% of the drift distance) and far from 
the wires (^^5% of the drift distance). The result¬ 
ing selection efficiency was 54% in BESS-Polar I 
and 56% in BESS-Polar IF 

5.4. Data Quality Cut Efficiencies 

The efficiencies of data quality cuts are deter¬ 
mined directly from flight data by the ratio of 
events after and before the cut, unlike the event 
selections discussed in Section lOl In BESS-Polar, 
these include track quality cuts and electronic haz¬ 
ard flag cuts. The rigidity of each particle is de¬ 
termined with good accuracy prior to these cuts. 
This is discussed in more detail in Section [5Tl 

To assure the accuracy of the rigidity measure¬ 
ment, the following track quality cuts were ap¬ 
plied: (1) reduced degree of freedom <5 
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for track fits in both rtf) and yz planes (where y is 
the vertical direction), (2) track fitting pathlength 
>500 mm, (3) ratio of fitted hits in the JET to ex¬ 
pected hits >0.6, (4) ratio of dropped hits in the 
JET and IDCs to expected hits <0.25, (5) 2 hits 
in each of the 2 IDCs all giving both x and z po¬ 
sitions, and (6) consistency between independent 
z positions measured by the tracking system and 
the TOE (determined using the time difference be¬ 
tween each end of the hit paddle). The efficiencies 
from (1) to (4) were ^ 98% for protons and helium 
nuclei in both flights. 

At 10 GeV/nucleon, the efficiencies for protons 
surviving (5) were (89.53 ± 0.08)% in BESS-Polar 
I and (90.37 ± 0.03)% in BESS-Polar II. The effi¬ 
ciencies from (5) for helium were (72.09±0.46)% in 
BESS-Polar I and (71.07 ± 0.18)% in BESS-Polar 
II. Helium efficiencies were almost 20% lower than 
for protons due to misidentification from delta- 
rays and saturation (overflow) of the FADCs caus¬ 
ing 5% of IDC hits to be cut. 

At 10 GeV/nucleon, the efficiencies for protons 
surviving (6) were (97.55 ±0.04)% for BESS-Polar 
I and (98.61 ± 0.01)% for BESS-Polar II. The effi¬ 
ciencies from (6) for helium were (98.26 ± 0.13)% 
for BESS-Polar I and (98.65 ± 0.04)% for BESS- 
Polar II. In BESS-Polar I, those detectors with 
only a single operating PMT did not provide re¬ 
dundant z positions and cut (6) could not be used. 

The TOE QDCs incorporated hazard flags that 
were set when the QDC was performing a func¬ 
tion that might give an incorrect measurement. 
Events displaying hazard flags were rejected. In 
BESS-Polar I, the QDGs were operated in a mode 
in which alternating capacitors were used to sam¬ 
ple background noise, which was then subtracted 
electronically from the measurements. This pro¬ 
cess proved to be inefficient at flight trigger rates. 
In BESS-Polar II, the efficiency was greatly im¬ 
proved because the QDGs were operated with a de¬ 
layed input and no background subtraction was re¬ 
quired. The efficiency of surviving this cut was es¬ 
timated as (52.00±0.09)% for protons and (52.23± 
0.35)% for helium nuclei at 10 GeV/nucleon in 
BESS-Polar I, and as (99.62 ± 0.01)% for protons 
and (99.62 ± 0.02)% for helium nuclei in BESS- 
Polar II. The trigger live time during which events 
were recorded was measured exactly by counting 1 
MHz clock pulses as 472881 seconds in BESS-Polar 
I and 1273381 seconds in BESS-Polar 11. Gombin- 


ing the trigger live time and the hazard flag cut 
efficiency give the actual measurement live time 
Tiive used in the flux calculations. 

5.5. Flight Data Estimated EfRciencies 

The final class of efficiencies that must be con¬ 
sidered are due to instrument operation or event 
analysis, but do not result from event selections. 
These are estimated from independent flight data 
selections. 

The event trigger efficiencies depend on the de¬ 
tection efficiencies from discrimination of signals 
in the UTOF and LTOF paddles through which 
a particle passed and the efficiency of the event 
trigger logic. These were estimated using parti¬ 
cles that triggered coincidences of the UTOF and 
MTOF, passed through all TOE layers (UTOF, 
MTOF, LTOF), and had i? > 10 GV. The trig¬ 
ger efficiencies were estimated from the probabil¬ 
ity that a normal event trigger was also gener¬ 
ated. These were (99.6 ± 0.1)% in BESS-Polar 1 
and (99.7 ± 0.1)% in BESS-Polar 11. 

The JET chamber has insensitive regions very 
close to the sense wires. When a particle passes 
through several such regions near the center of the 
JET, the track reconstruction algorithm cannot 
succeed. The corresponding track reconstruction 
efficiencies were estimated by scanning 5 x 10^ 
unbiased events with appropriate TOF hits and 
enough overall JET points. The efficiencies for 
both proton and helium were (99.0 ± 0.2)% in 
BESS-Polar 1 and (98.7 ± 0.2)% in BESS-Polar 
H. 

When a particle unrelated to a measured event 
(an “accidental”) passes through the JET or TOF 
within the resolving time, the single-track selec¬ 
tion discussed in Section [57^ eliminates the event. 
The influence of such accidental particles can be 
estimated by utilizing data recorded for random 
triggers that were issued independent of particle 
triggers at a rate of 10 Hz throughout the flight. 
Particles seen in these random trigger events were 
used to estimate the rates of accidental particles 
that might have caused real events to fail the 
single-track selections. The accidental efficiencies 
were 97.7% in BESS-Polar 1 and 96.2% in BESS- 
Polar 11. 
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Fig. 3.— Proton and helium bands in dE/dx vs 
R (top: UTOF, middle: LTOF, bottom: JET) 
from BESS-Polar 11. The superimposed lines show 
the selection criteria for protons and helium nu¬ 
clei. The inset graphs show the results of those 
selections above 30 GV. 

5.6. Particle Identification 

Protons and helium nuclei were identified by 
dE/dx vs R and /3“^ vs R techniques using the 
UTOF and LTOF. As examples, Figures |3] and 
m show the selection bands for protons and he¬ 
lium in BESS-Polar II. In BESS-Polar I, protons 
and helium nuclei were identified in the same man¬ 
ner. The efficiency of the TOF dE/dx vs R selec¬ 
tions were estimated by applying them to an inde¬ 
pendent sample of events selected using the JET 
dE/dx vs R. Overall, (95.84 ± 0.50)% of protons 
and (92.19 ± 1.50)% of helium nuclei were prop¬ 
erly identified in BESS-Polar I. In BESS-Polar II, 
(96.09 ±1.00)% of protons and (93.39 ± 1.50)% of 



Fig. 4.— BESS-Polar II proton (top) and helium 
(bottom) bands in /3~^ vs R after the selections 
shown in Figure O The superimposed lines show 
the selection criteria for protons and helium. The 
inset graphs show the results of those selections 
above 30 GV. 

helium nuclei were properly identified. The errors 
were estimated by comparing the TOF dE/dx ef¬ 
ficiencies in the flight data and in the simulated 
data. 

Since the /3“^ distribution is a Gaussian and the 
half-width of the selection band was set at 4 

(7, the efficiency is very close to unity for protons. 
In order to retain helium nuclei which might have 
been misidentified due to the contributions of ac¬ 
companying (5-rays, the width of the selection band 
was wider in the helium analysis at 8 a. With this 
adjustment the helium selection efficiency was also 
close to unity. 

5.7. Measurement Dependence on Rigid¬ 
ity 

The proton and helium measurements made by 
the BESS-Polar instruments span wide ranges of 
R. Their uniform performance over these wide 
ranges is demonstrated by the very small depen¬ 
dence on R of two figures-of-merit, the overall 
detection efficiency and the effective acceptance. 


























































































They are both nearly constant, especially at higher 
R, giving great confidence in the measured slopes 
of the interstellar proton and helium spectra. 

The track quality cuts required the positions 
hits in the tracking detector used in the analysis 
to fall near the fitted track. Analysis of these cuts 
at rigidities where multiple scattering can be ne¬ 
glected show that they do not depend on rigidity. 
Below 100 GV the tracking resolution before track¬ 
ing quality cuts is adequate to measure the rigidity 
to a small fraction of a rigidity bin and fluxes with 
and without quality cuts are essentially identical. 
There is no migration of events between rigidity 
bins due to the cuts. 

Figure [S] (top) shows the detection efficiencies, 
Edet, for protons and helium nuclei from the BESS- 
Polar II analysis, combining the track quality, 
track reconstruction, trigger, accidental, and par¬ 
ticle identification efficiencies. At the lowest R, 
some signals in the JET overflowed the FADCs due 
to high dE/dx and inconsistent z tracking between 
the JET and IDC slightly reduced detection effi¬ 
ciencies for helium. Above ^2 GV the efficiencies 
for both protons and helium were essentially con¬ 
stant. The detection efficiencies were stable with 
time and the efficiencies above 20 GV for 4-hour 
time intervals are estimated as (78.62 ± 0.93)% 
for protons and (59.35 ± 1.88)% for helium nu¬ 
clei in BESS-Polar I. In BESS-Polar II, they were 
(78.21 ± 0.60)% for protons and (59.41 ± 1.21)% 
for helium nuclei. 

To illustrate the uniformity of the BESS-Polar 
particle acceptance over the full range of R, we de¬ 
fine an effective acceptance as the product of the 
geometric acceptance (STl) and the selection effi¬ 
ciencies discussed in Section [5.31 This is shown in 
Figure[5] (bottom). The effective acceptances at 10 
GeV/nucleon in BESS-Polar I are 0.043 m^sr for 
protons and 0.039 m^sr for helium. In BESS-Polar 
II they are 0.058 m^sr for protons and 0.052 m^sr 
for helium. The systematic error in the effective 
acceptance primarily results from the systematic 
error in the non-interaction efficiency. 

5.8. Backgrounds and Measured Species 
5.8.1. Light Particles 
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Fig. 5.— BESS-Polar II detection efficiencies for 
protons and helium (top) and effective acceptances 
(bottom) as functions of R. 


duce light particles, including positrons, pions, 
and muons, that form a background to the pro¬ 
ton measurements. Below 1.6 GV the protons 
are clearly separated from the lighter particles, as 
shown in Figure IH and so are uncontaminated. 
However, above 1.6 GV the light particles can no 
longer be identified. To estimate the magnitude 
of the resulting background, the ratio of light par¬ 
ticles to protons was measured as a function of 
rigidity in the range where both can be clearly 
identified. The light particle to proton ratio was 
2.1% at 0.8 GV, rapidly dropping to 0.8% at 1.6 
GV, and could be fit well by a power law in R with 
an index of —1.6, indicating that contamination is 
negligible above 1.6 GV. The measurements were 
reproduced well by solving simultaneous trans¬ 
port equations using as input BESS-Polar II pro¬ 
tons and helium nuclei and electrons and positrons 
from PAMELA measurements. Because both the 
measurements and calculations showed that con¬ 
tamination of the measured protons by light par¬ 
ticles is negligible, no corrections were made. 


Interactions of cosmic rays with the atmosphere 
above BESS-Polar and with the balloon-flight sus¬ 
pension components, especially the swivel, pro- 
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5.8.2. Deuterons and Tritons 

The BESS-Polar spectrometers distinguish 
deuterons up to 3 GV and our measurements 
show that the d/p ratio decreases with increasing 
rigidity, falling to 2.5% at 3 GV in BESS-Polar 
1. Because of its greater atmospheric overburden, 
the d/p ratio at 3 GV is 3.0% in BESS-Polar II. 
Above 3 GV, the deuterons could not be resolved. 
Similarly, tritons can be resolved up to 5 GV, 
where the measured t/p ratio was ~ 0.5%. Above 
5 GV tritons could not be resolved. Thus, above 3 
GV the reported protons have a small admixture 
of deuterons and above 5 GV an even smaller ad¬ 
mixture of tritons. We conservatively include '^3% 
uncertainty for deuteron contamination that may 
have increased the measured proton flux. Tritons 
were not included in the uncertainty. 

5.8.3. ^He 

In both flights, Z=2 nuclei were clearly identi¬ 
fied by dE/dx vs R and vs R techniques using 
the UTOF and LTOF hodoscopes and were con¬ 
firmed using dE/dx vs R from the JET. However, 
these nuclei included both ^He and ^He. Follow¬ 
ing the practice used in previous BESS reports of 
He fluxes, all Z=2 particles were included in the 
He spectrum. 


5.9. Flux Calculations and Corrections 

5.9.1. Top of Instrument Flux 

Nuclei that have enough energy to penetrate to 
the LTOF and trigger the instrument have nearly 
constant energy throughout the tracker. Thus, the 
energy that the particle would have had when en¬ 
tering the instrument (top-of-instrument or TOI) 
is calculated by integrating its differential energy 
losses dE/dx for all materials above the tracker 
using a detailed model of the instrument compo¬ 
nents encountered along its trajectory. The mea¬ 
sured differential flux of a given particle species at 
TOI (d>TOi)j binned in energy with width dE, can 
then be determined by: 


d’TOl(ATOl) 




^det * ^MC ‘ SH • Tji, 


( 2 ) 


where Np is number of observed particles of that 
species. 


5.9.2. Spectral Deconvolution 

For steeply falling particle spectra, observed 
spectral shapes can be deformed because of the 
R dependence of the magnetic-rigidity spectrome¬ 
ter R resolution, or because of the energy depen¬ 
dence of the energy resolution of a calorimeter, 
extensive air shower detector, etc.. The details 
of this deformation depend on the absolute R or 
absolute energy dependence of the measurement 
resolution. Response functions of different mea¬ 
surement techniques differ with energy and some¬ 
times (e.g. calorimeters) with particle species. For 
magnetic-rigidity spectrometers like BESS-Polar, 
the resolution deteriorates as R approaches the 
MDR. The more abundant lower R particles are 
misidentified as higher i?, and, by number con¬ 
servation, fluxes are correspondingly depleted at 
slightly lower R, giving an apparent dip in the 
spectrum and then a recovery to a harder spectral 
index than the underlying spectrum. Depending 
on the underlying spectral index, small R errors of 
this sort can result in a considerable change in the 
observed spectral shape and power law index (or 
indices). To compensate, reported high R spectra 
from magnetic rigidity spectrometers are usually 
the result of a “deconvolution” procedure to repro¬ 
duce the most probable real spectrum that would 
have resulted in the observed spectrum, taking 
into account the finite R resolution. For the re¬ 
sults reported here, the deconvolution was accom¬ 
plished using a statistical procedure based on iter¬ 
atively applying Bayes’ theorem. In this method, 
the finite resolution effect is described by a con¬ 
volution matrix in deflection space, which was es¬ 
timated from a simulation using the appropriate 
(i?“^) resolution. As noted above, for both BESS- 
Polar I and BESS-Polar H, (j{R-^) = 5.0 TV"!. 
An initial spectrum must be chosen as the starting 
point for the deconvolution. Based on results from 
previous measurements in the BESS-Polar rigid¬ 
ity range, single power laws with spectral indices 
of —2.80 for protons and —2.70 for helium nuclei 
were used. The systematic error of the deconvo¬ 
lution procedure was estimated to be about 1% 
at 160 GV from the deviation due to input index 
errors (±0.05) and (j{R~^) error (±0.3 TV“^). 


10 



5.9.3. Atmospheric Loss and Production Correc¬ 
tions 


In order to obtain the absolute flux of pri¬ 
mary protons and helium nuclei at TOA, correc¬ 
tions were applied for the survival probability in 
the residual atmosphere ( 77 ) and estimated atmo¬ 
spheric secondary production ratio (i?air)- Both 
were estim ated by solving sim ultaneous transport 
equations ( Papini et alJll996r) . The primary spec¬ 
trum at TOA ($toa) was determined in an iter¬ 
ative procedure so that the estimated spectrum 
at TOI ($TOi) agreed with the observed spec¬ 
trum. For BESS-Polar I, with residual atmo¬ 
sphere of 4.33 g/cm^, the survival probability at 
10 GeV/nucleon was estimated as (94.42 ± 0.56)% 
for protons and (90.77 ± 1.2)% for helium. For 
BESS-Polar II, with residual atmosphere of 5.81 
g/cm^, the survival probability was estimated as 
(92.47 ± 0.74)% for protons and (87.63 ± 1.5)% 
for helium. The proton secondary-to-primary ra¬ 
tio at 10 GeV was calculated as 0.0174 ±0.0020 in 
BESS-Polar 1 and 0.0227 ± 0.0026 in BESS-Polar 
11. Atmospheric secondary helium nuclei above 
1 GeV/nucleon are dominated by the fragments 
of heavier cosmic-ray nuclei. The secondary-to- 
primary helium ratio at 10 GeV/nucleon was cal¬ 
culated as 0.0209 ± 0.0043 in BESS-Polar 1 and 
0.0271 ± 0.0054 in BESS-Polar 11. The errors in 
atmospheric correction were estimated by multi¬ 
plying the factors which represent the uncertainty 
of residual air depth as 10%, cross section of pri¬ 
mary protons as 5%, and the cross section of pri¬ 
mary heavier cosmic ray nuclei as 20%. 


5 . 9 . 4 . Top of Atmosphere Flux 

Utilizing these atmospheric loss and production 
calculations, the measured differential flux of a 
given particle species extrapolated at the top of 
the atmosphere (4>toa) and binned in energy with 
width dE is expressed as: 


corresponding Eatm- 

6. Systematic Uncertainties 

The overall systematic uncertainties for protons 
and helium share nine components; particle iden¬ 
tification, track reconstruction, trigger, geomet¬ 
ric acceptance, non-interacted single track selec¬ 
tion, atmospheric production and loss correction, 
spectral deconvolution, magnetic field details, and 
tracking system alignment. For protons, deuteron 
contamination may also contribute by increasing 
the measured proton flux. Details of those factors 
are discussed above and their relative values as 
functions of particle energy are illustrated in Fig¬ 
ure [S] which shows the systematic and statistical 
uncertainties in BESS-Polar II proton measure¬ 
ments. The total systematic uncertainty in the 
proton measurements is asymmetric and is shown 
both with and without deuteron contamination. 
The upper (additive) total systematic uncertainty 
in the measured proton spectrum does not include 
deuteron contamination. The lower (subtractive) 
total systematic uncertainty in the measured pro¬ 
ton spectrum includes deuteron contamination. 
The proton total systematic uncertainties at 10 
GeV are 1.9% (upper), 3.4% (lower) in BESS- 
Polar 1 and 2.1% (upper), 3.7% (lower) in BESS- 
Polar 11. The difference arises from the greater at¬ 
mospheric overburden in BESS-Polar 11. The he¬ 
lium total systematic uncertainties are symmetric 
above and below the measured spectral points. At 
10 GeV/nucleon they are 3.0% in BESS-Polar 1, 
and 3.2% in BESS-Polar 11. However, the uncer¬ 
tainty of the non-interacted single track selection 
for helium is worse at high energy in BESS-Polar 
I because it is difficult to simulate the effect of d- 
rays and interactions for TOF paddles with only 
one PMT. Thus, at 80 GeV/nucleon the helium 
total systematic uncertainties are 5.3% in BESS- 
Polar 1, and 5.0% in BESS-Polar 11. 


$toa(Atoa) = — ij , - F -i 

9{TjATm) + Ra_ir{EATM) 

where Eatm is a shorthand for the energy used 
in each step of the iterative process described 
above in which the energy of the particle pass¬ 
ing through the atmosphere is calculated at 0.1 
g/cm? intervals. The loss and production calcula¬ 
tions are made over the same intervals using the 


7. Measured Spectra 

We have measured the absolute TOA fluxes 
of protons 0.2-160 GeV and helium nuclei 0.15- 
80 GeV/nucleon using both BESS-Polar 1 and 
BESS-Polar 11, where the lower energy was set by 
the requirement for particles to penetrate to the 
LTOF and upper energy was chosen to give flat 
detection efficiency. The corresponding ranges are 
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Fig. 6.— Individual and total systematic uncer¬ 
tainties and the statistical uncertainty for BESS- 
Polar II proton measurements. The total system¬ 
atic uncertainty for protons is asymmetric. The 
lower (subtractive) uncertainty includes deuteron 
contamination but the upper (additive) does not. 


0.6-160 GV for protons and 1.1-160 GV for helium. 
The BESS-Polar measurements are summarized in 
Tables [T] and [5] The overall uncertainties, includ¬ 
ing both statistical and systematic errors, are less 
than ±10% for BESS-Polar I protons, ±6% for 
BESS-Polar II protons, ±7% for BESS-Polar I he¬ 
lium, and ±5% for BESS-Polar II helium. The 



took place at very low and nearly constant geo¬ 
magnetic cutoff. Their measured spectra do not 
have the systematic uncertainties at lower ener¬ 
gies that result from the rigidity-dependent live 
time required by orbiting instruments such as 
PAMELA (70° inclination), and AMS-02 (51.6° 
inclination) to account for measurements in vary¬ 
ing geomagne tic cutoff. GREA M proton and he¬ 
lium spectra ( Ahn et al. 2011)11 are not shown in 
Figure [7] because most points are at higher ener- 


7.1. High Energy Interstellar Spectra 

At energies above 30 GeV for protons and 15 
GeV/nucleon for helium, the effects of solar mod¬ 
ulation (discussed in Section 17.2() are effectively 
negligible and the proton spectra measured BESS- 
Polar I and BESS-Polar II are essentially identi¬ 
cal, differing by less than 1% at 160 GeV despite 
having been obtained during quite different pe¬ 
riods of solar activity (see further discussion be¬ 
low). The same is true of the helium spectra 
above ^15 GeV/nucleon where the BESS-Polar I 
and BESS-Polar II results differ by less than 1% 
at 80 GeV/nucleon. Thus, the interstellar spec¬ 
tra can be assumed to be the same as the mea¬ 
sured TOA spectra over 30-160 GeV for protons 
and 15-80 GeV/nucleon for helium. Over these 
ranges, the measured spectra appear to follow sin¬ 
gle power laws. Thus, a given spectral F can be 
parameterized by a power law in kinetic energy, 

Ek, as F(Ek) = ^Ek^k ^ where ^Ek is the nor¬ 
malization constant. For the results reported here, 
systematic errors related to spectral index (decon¬ 
volution and alignment) were taken into account, 
and a fitting range of 30-160 GeV was used for 
protons and 15-80 GeV/nucleon for helium. The 
best fit values and uncertainties for protons (q^*") 
and helium nuclei (y^e) for the two flights are: 

7f(kiari) = - 2-715 ± 0.006(stat) ± O.OlO(syst) 
7f(Poiarii) = - 2-736 ± 0.002(stat) ± O.OlO(syst) 
7He(Poiari) = “ 2-606 ± 0.017(stat) ± 0.012(syst) 
7HeVoiarii) = “ 2.627 ± 0.006(stat) ± 0 . 012 (syst) 

For both protons and helium, fits can be made 
over 30-160 GV to a single power law in i?, pa¬ 
rameterized as F(R) = . The best fit 

values and uncertainties over a fitting range of 30- 
160 GV for both protons ( 7 ^) helium (y/ig) for the 
two flights are: 

7^(Poiari) = “ 2.763 ± 0.006(stat) ± O.OlO(syst) 
7^(Poiarii) = - 2-784 ± 0.002(stat) ± O.OlO(syst) 
7He(Poiari) = “ 2-690 ± 0.016(stat) ± O.Oll(syst) 
7He(Poiarii) = “ 2.709 ± 0.006(stat) ± O.Oll(syst) 
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Fig. 7.— Absolute differential energy spectra of primary protons and helium nuclei at TQA multiplie d 
by me asured by BESS-Pol ar I and B ESS-Polar II Spec tra reported by BESS98 dSanuki et alJ 2Q0Q|) . 
BESS-Te V (iHaino et all l2004a[l . ATIC-2 dPanov et al.ll2007l) . PAMELA (|Adriani et alJl201ll l2014l) . and 
AMS-02 ( Aguilar et alJ2015ai b) are shown for comparison. Differences in the BESS-Polar I and BESS-Polar 
II spectra at low energies are due to solar modulation. 


At high-ener gy, PAMELA (lAdria.ni et al. 2011 


2014 ( 1 . AMS-02 (lAgnilar et al.H2015alllJi andBESS- 


Polar proton and helium spectra agree within one 
aE- We note, however, that the BESS-Polar spec¬ 
tra do not extend to high enough energy to test the 
spectral harde ning re ported by PAMELA, A TIC- 
2 ()Panov et al.ll2007ll . CREAM ()Ahn et al.ll201^ 
and AMS-02. 


7.2. Low Energies and Solar Modulation 

GCR entering the heliosphere are scattered by 
irregularities in the heliospheric magnetic held 


(HMF) and undergo convection and adiabatic de¬ 
celeration in the expanding solar wind. The re¬ 
sulting modihcation of the GCR energy spectra, 
is known as “solar modulation”. Rehecting varia¬ 
tions in solar activity and HMF polarity, large dif¬ 
ferences at low energies have been observed in the 
proton and helium spectra measured by the BESS 
collaboration since 1993. The solar cycle includes 
both a sunspot activity period of ^^11 years and a 
^22 year magnetic cycle with alternating positive 
(A>0) and negative (A<0) phases. A<0 polarity 
cycles are dehned as the periods when the HMF 
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Fig. 8.— Variation of Bartol neutron monitor 
counts and sunspot number from 1991 to 2015. 
The BESS, ATIC-2 and BESS-Polar flights are 
marked as circles at the corresponding sunspot 
number. The periods for which PAMELA and 
AMS-02 proton and helium spectra have been re¬ 
ported are shown as bars at the approximate aver¬ 
age sunspot number. The BESS-Polar II flight was 
carried out very near the deepest solar minimum. 


is directed towards the Sun in its northern hemi¬ 
sphere. HMF polarity reverses at a time when 
solar activity is maximum, and the global mag¬ 
netic field profile also reverses throughout the he¬ 
liosphere with time lag due to propagation of the 
fields. Positive and negative particles traversing 
the HMF drift in opposite directions, taking dif¬ 
ferent routes to arrive at the Earth and giving rise 
to “flat” and “peaked” periods in neutron monitor 
data around “solar minimum”. Thus, the details 
of the effect of the solar wind and its entrained 
magnetic fields on the incoming GCR have to be 
taken into account in deriving interstellar spec¬ 
tra at lower energies. At energies below 30 GeV 
for protons and 15 GeV/nucleon for helium (dis¬ 
cussed above) we report only the measured TOA 
fluxes. However, we note the rich information on 
solar modulation that can be derived from BESS 
and BESS-Polar antiproton, proton, helium, and 
light isotope spectra obtained over more than a 
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Fig. 9.— Absolute differential energy spectra of 
primary protons measured by BESS-Polar I and 
BESS-Polar H together with earlier BESS inea- 
sureme nts and the PAME LA spectra compiled 
yearly ( Adrian! et al.l [201^ . The effect of solar 
modulation is evident at lower energies. Low en¬ 
ergy measurements at similar solar modulation, 
e.g. BESS97, BESS-Polar H, and PAMELA07 are 
consistent within reported errors. The AMS-02 re¬ 
sults in Figure [7] are not reported by time period. 


full solar cycle, including a reversal of solar mag¬ 
netic polarity. To illustrate the range of solar 
activity spanned by BESS and BESS-Polar, Eig- 
ure [5] shows the Bartol neutr on monitor c ount- 
ing rate (blue dashed curve) (iBiebed l20I4ll and 
the n umber of sunspots (orange curve) (iHathaw^ 
2 OI 5 II together the data of BESS flights (red cir¬ 
cles and blue stars). The ATIC-2 (light blue cir¬ 
cle) , PAMELA (gray bar) and AMS-02 (green bar) 
flights are also indicated. 

Differences in the spectra measured by BESS- 
Polar I and BESS-Polar H from solar activity be¬ 
low ~10 GeV for protons in Figure [9] and below 
^5 GeV/nucleon for helium in Figure [TUI are high- 
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Fig. 10.— Absolute differential energy spectra 
of helium nuclei measured by BESS-Polar I and 
BESS-Polar II together with earlier BESS mea¬ 
surements. BESS97 and BESS-Polar II measure¬ 
ments at similar solar modulation are consistent 
within reported errors except at the lowest ener¬ 
gies. 


lighted and compared to measurements by ear¬ 
lier BESS instruments and to PAMELA measure¬ 
ments for four years before and during the most 
recent solar minimum. Varying solar modulation 
strongly affected proton and helium spectra mea¬ 
sured by BESS (BESS 97, 98, 99) approaching 
the HMF polarity reversal from positive to neg¬ 
ative in 2000, prior to the BESS 00 flight. Proton 
and helium measurements following the HMF po¬ 
larity reversal (BESS 00, BESS-TeV, BESS-Polar 
I, BESS-Polar II) show much smaller variations 
with solar activity. This is consistent with the 
effect expected from charge-sign dependent drift 
(iBieber et al. I999I : Asaoka et al. 2002t IPotgieter 


20141 1. The proton spectra observed by BESS at 


solar minimum in 1997 and and by BESS-Polar II 
and PAMELA in 2007, near the next solar mini¬ 


Fig. 11.— Proton/helium flux ratios measured 
by BESS-Polar I and BESS-Polar II (with pro¬ 
tons rebinned to mat ch helium). The ratio re - 
ported by PAMELA ( Adriani et al.ll^01lL 12014 1 
and the AMS-02 ratio ( Aguilar et al.l ~15b[l are 
also shown. Above 10 GV the two BESS-Polar 
ratios and AMS-02 agree well within error bars, 
while the PAMELA ratio is slightly higher and 
flatter with increasing R. 


mum are very similar. As would be expected, the 
low-energy PAMELA proton spectrum from 2006, 
prior to solar minimum falls below BESS-Polar II 
or PAMELA in 2007, and the PAMELA results 
from 2008 and 2009 during the deepest solar min¬ 
imum are higher. The BESS and BESS-Polar mea¬ 
surements of antiprotons and protons, which differ 
only in charge-sign and interstellar spectral shape, 
provide an excellent test c ase for charge-sign de¬ 
pendent solar modulation ( Asaoka et al.l 2002 1. A 
full discussion solar modulation and of the impli¬ 
cations of the BESS results is beyond the scope of 
the present paper. 
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7.3. Proton-to-Helium Ratio 


Comparing the spectra of protons and helium 
gives important information on the sources and 
acceleration of the light cosmic-rays and a test of 
cosmic-ray origin and propagation models. The 
measured power-law indices of protons and he¬ 
lium, given above, are slightly different. This is 
most easily shown in R space because our measure¬ 
ments for both species reach 160 GV. The mea¬ 
sured proton/helium flux ratios in R from BESS- 
Polar I and BESS-Polar II, with protons re-binned 
to match helium, are s hown in Figure TT]an d com - 
pared to PAMELA ( Adriani et al. 2 0IlL 20141 
data and with AMS-02 results (|Aguilar et al 
2015bll . The BESS-Polar proton/helium flux ra¬ 


tios are summarized in Table [3l A simple model 
of solar modulation, appropriate when comparing 
two species of the same charge-sign, shows that the 
effect of solar modulation on this ratio should be 
negligible above 20 GV. This is borne out by the 
measured ratios. As expected, p/He above 20 GV 
is well represented by a single power-law with an 
index which is just the difference of the indices of 
the protons and helium as = 7^ —7He)- This 
was also determined directly using a power-law fit 
to p/He above 20 GV, as shown in Figure [TT] 
The fit values obtained were A^n-fu = —0.081 ± 
O.OlO(stat) ± 0.020(syst) for BESS-Polar I and 
A^R-fit = —0.086 ± 0.004(stat) ± 0.022(syst) for 
BESS Polar II. These compare well to the calcu¬ 
lated values of A^r-cuIc = —0.073 ± 0.017(stat) ± 
0.015(syst) for BESS-Polar I and A^a-cau = 
-0.075 ± 0.006(stat) ± 0.015(syst) for BESS-Polar 
II. Above 10 GV the two BESS-Polar ratios and 
AMS-02 agree well within error bars, while the 
PAMELA ratio is slightly higher and flatter with 
increasing R. 


8. Conclusion 

With the BESS-Polar I long-duration balloon 
flight in 2004 and the BESS-Polar II flight in 
2007/2008, we have measured the energy spectra 
at the top of the atmosphere of protons in the 
range 0.2-160 GeV and helium nuclei in the range 
0.15-80 GeV/nucleon, corresponding to magnetic- 
rigidity ranges of 0.6-160 GV for protons and 1.1- 
160 GV for helium. At energies above the effects 
of solar modulation, the spectra of either species 
measured in the two flights differ by less than 1% 


at the highest measured energy (rigidity). Within 
the measured ranges, the high energy (rigidity) 
spectra of both species can be represented very 
well by single power laws. We have also deter¬ 
mined the measured proton to helium ratio and 
show that at high energies it is also fit by a single 
power law. 
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Table 1 

Proton flux at the top of atmosphere 




BESS-Polar I 

BESS-Polar II 

Energy Range 

Ek 

Flux ± AFsta ± AFsys 

Flux ± AFsta ± AFsys 

(GeV) 

(GeV) 

(m= sr s GeV)“^ 

(m^ sr s GeV)“^ 


0 . 200 - 0.211 

0.205 

4.825 

± 

0.032 

± 

0.446 

X 

10 ^ 

1.527 

± 

0.002 

± 

0.089 

X 

lo ' 

0 . 211 - 0.224 

0.218 

5.110 

± 

0.031 

± 

0.422 

X 

10 ^ 

1.567 

± 

0.002 

± 

0.086 

X 

10 ' 

0 . 224 - 0.237 

0.231 

5.353 

± 

0.030 

± 

0.423 

X 

10 ^ 

1.603 

± 

0.002 

± 

0.086 

X 

10 ' 

0 . 237 - 0.251 

0.244 

5.565 

± 

0.029 

± 

0.397 

X 

10 ^ 

1.633 

± 

0.002 

± 

0.082 

X 

10 - 

0 . 251 - 0.266 

0.259 

5.780 

± 

0.028 

± 

0.374 

X 

10 ^ 

1.662 

± 

0.002 

± 

0.079 

X 

10 - 

0 . 266 - 0.282 

0.274 

6.000 

± 

0.027 

± 

0.354 

X 

10 ^ 

1.684 

± 

0.001 

± 

0.076 

X 

10 - 

0 . 282 - 0.298 

0.290 

6.263 

± 

0.026 

± 

0.358 

X 

10 ^ 

1.702 

± 

0.001 

± 

0.075 

X 

10 ' 

0 . 298 - 0.316 

0.307 

6.399 

± 

0.025 

± 

0.336 

X 

10 ^ 

1.713 

± 

0.001 

± 

0.072 

X 

10 ' 

0 . 316 - 0.335 

0.326 

6.649 

± 

0.024 

± 

0.322 

X 

10 ^ 

1.722 

± 

0.001 

± 

0.069 

X 

10 ' 

0 . 335 - 0.355 

0.345 

6.787 

± 

0.024 

± 

0.304 

X 

10 ^ 

1.723 

± 

0.001 

± 

0.066 

X 

10 ' 

0 . 355 - 0.376 

0.365 

6.910 

± 

0.023 

± 

0.303 

X 

10 ^ 

1.718 

± 

0.001 

± 

0.065 

X 

10 ' 

0 . 376 - 0.398 

0.387 

7.061 

± 

0.022 

± 

0.288 

X 

10 ^ 

1.708 

± 

0.001 

± 

0.062 

X 

10 ' 

0 . 398 - 0.422 

0.410 

7.191 

± 

0.021 

± 

0.275 

X 

10 ^ 

1.693 

± 

0.001 

± 

0.059 

X 

10 ' 

0 . 422 - 0.447 

0.434 

7.319 

± 

0.021 

± 

0.262 

X 

10 ^ 

1.673 

± 

0.001 

± 

0.056 

X 

10 ' 

0 . 447 - 0.473 

0.460 

7.326 

± 

0.020 

± 

0.259 

X 

10 ^ 

1.645 

± 

0.001 

± 

0.055 

X 

10 ' 

0 . 473 - 0.501 

0.487 

7.373 

± 

0.019 

± 

0.246 

X 

10 ^ 

1.615 

± 

0.001 

± 

0.052 

X 

10 ' 

0 . 501 - 0.531 

0.516 

7.386 

± 

0.019 

± 

0.234 

X 

10 ^ 

1.582 

± 

0.001 

± 

0.050 

X 

10 ' 

0 . 531 - 0.562 

0.546 

7.391 

± 

0.018 

± 

0.222 

X 

10 ^ 

1.542 

± 

0.001 

± 

0.047 

X 

10 ' 

0 . 562 - 0.596 

0.579 

7.351 

± 

0.017 

± 

0.219 

X 

10 ^ 

1.501 

± 

0.001 

± 

0.045 

X 

10 ' 

0 . 596 - 0.631 

0.613 

7.295 

± 

0.017 

± 

0.208 

X 

10 ^ 

1.456 

± 

0.001 

± 

0.043 

X 

10 ' 

0 . 631 - 0.668 

0.649 

7.189 

± 

0.016 

± 

0.196 

X 

10 ^ 

1.411 

± 

0.001 

± 

0.040 

X 

10 ' 

0 . 668 - 0.708 

0.688 

7.116 

± 

0.015 

± 

0.187 

X 

10 ^ 

1.360 

± 

0.001 

± 

0.038 

X 

10 ' 

0 . 708 - 0.750 

0.729 

6.994 

± 

0.015 

± 

0.182 

X 

10 ^ 

1.308 

± 

0.001 

± 

0.036 

X 

10 ' 

0 . 750 - 0.794 

0.772 

6.862 

± 

0.014 

± 

0.172 

X 

10 ^ 

1.256 

± 

0.001 

± 

0.034 

X 

10 ' 

0 . 794 - 0.841 

0.818 

6.746 

± 

0.013 

± 

0.163 

X 

10 ^ 

1.207 

± 

0.001 

± 

0.032 

X 

10 ' 

0 . 841 - 0.891 

0.866 

6.571 

± 

0.013 

± 

0.154 

X 

10 ^ 

1.153 

± 

0.001 

± 

0.029 

X 

10 ' 

0 . 891 - 0.944 

0.917 

6.388 

± 

0.012 

± 

0.149 

X 

10 ^ 

1.099 

± 

0.001 

± 

0.028 

X 

10 ' 

0 . 944 - 1.000 

0.972 

6.218 

± 

0.012 

± 

0.142 

X 

10 ^ 

1.041 

± 

0.001 

± 

0.026 

X 

10 ' 

1 . 000 - 1.059 

1.029 

6.011 

± 

0.011 

± 

0.134 

X 

10 ^ 

9.875 

± 

0.005 

± 

0.243 

X 

10 - 

1 . 059 - 1.122 

1.090 

5.760 

± 

0.011 

± 

0.125 

X 

10 ^ 

9.327 

± 

0.005 

± 

0.225 

X 

lO ' 

1 . 122 - 1.189 

1.155 

5.535 

± 

0.010 

± 

0.120 

X 

10 ^ 

8.798 

± 

0.005 

± 

0.212 

X 

lO ' 

1 . 189 - 1.259 

1.223 

5.336 

± 

0.010 

± 

0.114 

X 

10 ^ 

8.276 

± 

0.004 

± 

0.197 

X 

lO ' 

1 . 259 - 1.334 

1.296 

5.078 

± 

0.009 

± 

0.107 

X 

10 ^ 

7.766 

± 

0.004 

± 

0.182 

X 

lO ' 

1 . 334 - 1.413 

1.373 

4.869 

± 

0.009 

± 

0.101 

X 

10 ^ 

7.270 

± 

0.004 

± 

0.168 

X 

10 - 

1 . 413 - 1.496 

1.454 

4.607 

± 

0.008 

± 

0.095 

X 

10 ^ 

6.794 

± 

0.004 

± 

0.157 

X 

10 - 

1 . 496 - 1.585 

1.540 

4.378 

± 

0.008 

± 

0.089 

X 

10 ^ 

6.333 

± 

0.003 

± 

0.145 

X 

10 - 

1 . 585 - 1.679 

1.631 

4.143 

± 

0.007 

± 

0.083 

X 

10 ^ 

5.896 

± 

0.003 

± 

0.133 

X 

10 - 

1 . 679 - 1.778 

1.728 

3.905 

± 

0.007 

± 

0.078 

X 

10 ^ 

5.473 

± 

0.003 

± 

0.123 

X 

10 - 

1 . 778 - 1.884 

1.830 

3.667 

± 

0.006 

± 

0.073 

X 

10 ^ 

5.065 

± 

0.003 

± 

0.113 

X 

10 - 

1 . 884 - 1.995 

1.939 

3.447 

± 

0.006 

± 

0.068 

X 

10 ^ 

4.673 

± 

0.003 

± 

0.104 

X 

10 - 

1 . 995 - 2.113 

2.054 

3.214 

± 

0.006 

± 

0.063 

X 

10 ^ 

4.307 

± 

0.002 

± 

0.095 

X 

10 - 

2 . 113 - 2.239 

2.175 

3.003 

± 

0.005 

± 

0.058 

X 

10 ^ 

3.956 

± 

0.002 

± 

0.086 

X 

10 - 

2 . 239 - 2.371 

2.304 

2.804 

± 

0.005 

± 

0.054 

X 

10 ^ 

3.628 

± 

0.002 

± 

0.079 

X 

10 - 

2 . 371 - 2.512 

2.441 

2.607 

± 

0.005 

± 

0.050 

X 

10 ^ 

3.318 

± 

0.002 

± 

0.072 

X 

10 - 

2 . 512 - 2.661 

2.585 

2.424 

± 

0.004 

± 

0.046 

X 

10 ^ 

3.041 

± 

0.002 

± 

0.066 

X 

lo ' 

2 . 661 - 2.818 

2.739 

2.242 

± 

0.004 

+ 0.043 
-0.043 

X 

10 ^ 

2.777 

± 

0.002 

+0.060 

-0.060 

X 

10 - 

2 . 818 - 2.985 

2.901 

2.059 

± 

0.004 

+ 0.039 
-0.039 

X 

10 ^ 

2.531 

± 

0.002 

+0.054 

-0.054 

X 

10 - 

2 . 985 - 3.162 

3.073 

1.906 

± 

0.004 

+ 0.036 
-0.037 

X 

10 ^ 

2.305 

± 

0.001 

+0.049 

-0.051 

X 

10 - 

3 . 162 - 3.415 

3.286 

1.709 

± 

0.003 

+ 0.032 
-0.034 

X 

10 ^ 

2.061 

± 

0.001 

+0.044 

-0.047 

X 

10 - 

3 . 415 - 3.687 

3.548 

1.513 

± 

0.003 

+ 0.028 
-0.033 

X 

10 ^ 

1.801 

± 

0.001 

+0.038 

-0.046 

X 

10 - 

3 . 687 - 3.981 

3.831 

1.330 

± 

0.002 

+ 0.025 
-0.032 

X 

10 ^ 

1.566 

± 

0.001 

+0.033 

-0.045 

X 

10 - 

3 . 981 - 4.299 

4.137 

1.164 

± 

0.002 

+ 0.022 
-0.029 

X 

10 ^ 

1.356 

± 

0.001 

+0.029 

-0.040 

X 

10 - 

4 . 299 - 4.642 

4.467 

1.015 

± 

0.002 

+ 0.019 
-0.026 

X 

10 ^ 

1.171 

± 

0.001 

+0.025 

-0.035 

X 

10 - 

4 . 642 - 5.012 

4.823 

8.816 

± 

0.017 

+ 0.163 
-0.240 

X 

10 ^ 

1.006 

± 

0.001 

+0.021 

-0.032 

X 

10 - 

5 . 012 - 5.412 

5.208 

7.612 

± 

0.015 

+ 0.141 
-0.221 

X 

10 ^ 

8.621 

± 

0.006 

+0.179 

-0.285 

X 

10 - 

5 . 412 - 5.843 

5.623 

6.557 

± 

0.013 

+ 0.121 
-0.198 

X 

10 ^ 

7.364 

± 

0.005 

+0.153 

-0.248 

X 

10 - 


18 






Table 1 —Continued 


BESS-Polar I BESS-Polar II 


Energy Range Ek Flux ± AF^ta ± AFgys Flux ± AF^ta ± AFgys 


( GeV ) 

( GeV ) 


( m ^ sr s 

GeV )" 

1 



( m ^ sr s 

GeV )" 

1 


5 . 843 - 6.310 

6.072 

5.608 

± 

0.012 

+ 0.103 
-0.176 

X 

10 ^ 

6.257 

± 

0.005 

+0.130 

-0.215 

X 

10 ^ 

6 . 310 - 6.813 

6.556 

4.790 

± 

0.010 

+ 0.088 
-0.151 

X 

10 ^ 

5.305 

± 

0.004 

+0.110 

-0.184 

X 

10 ^ 

6 . 813 - 7.356 

7.079 

4.083 

± 

0.009 

+ 0.075 
-0.130 

X 

10 ^ 

4.481 

± 

0.004 

+0.093 

-0.156 

X 

10 ^ 

7 . 356 - 7.943 

7.644 

3.449 

± 

0.008 

+ 0.063 
-0.112 

X 

10 ^ 

3.771 

± 

0.003 

+0.078 

-0.132 

X 

10 ^ 

7 . 943 - 8.577 

8.254 

2.930 

± 

0.007 

+ 0.054 
-0.096 

X 

10 ^ 

3.169 

± 

0.003 

+0.065 

-0.111 

X 

10 ^ 

8 . 577 - 9.261 

8.912 

2.457 

± 

0.006 

+ 0.045 
-0.081 

X 

10 ^ 

2.644 

± 

0.003 

+0.055 

-0.093 

X 

10 ^ 

9 . 261 - 10.00 

9.623 

2.050 

± 

0.006 

+ 0.037 
-0.067 

X 

10 ^ 

2.211 

± 

0.002 

+0.045 

-0.077 

X 

10 ^ 

10 . 00 - 11.01 

10.49 

1.682 

± 

0.004 

+ 0.031 
-0.055 

X 

10 ^ 

1.797 

± 

0.002 

+0.037 

-0.063 

X 

10 ^ 

11 . 01 - 12.12 

11.55 

1.336 

± 

0.004 

+ 0.024 
-0.044 

X 

10 ^ 

1.419 

± 

0.001 

+0.029 

-0.050 

X 

10 ^ 

12 . 12 - 13.34 

12.71 

1.059 

± 

0.003 

+ 0.019 
-0.035 

X 

10 ^ 

1.123 

± 

0.001 

+0.023 

-0.039 

X 

10 ^ 

13 . 34 - 14.68 

13.99 

8.393 

± 

0.027 

+ 0.154 
-0.277 



8.802 

± 

0.010 

+0.182 

-0.309 



14 . 68 - 16.16 

15.40 

6.563 

± 

0.023 

+ 0.121 
-0.216 



6.910 

± 

0.009 

+0.143 

-0.243 



16 . 16 - 17.78 

16.95 

5.156 

± 

0.019 

+ 0.095 
-0.170 



5.393 

± 

0.007 

+0.112 

-0.190 



17 . 78 - 19.57 

18.65 

4.059 

± 

0.016 

+ 0.076 



4.201 

± 

0.006 

+0.088 



19 . 57 - 21.54 

20.53 

3.147 

± 

0.014 

+ 0.059 



3.271 

± 

0.005 

+0.069 



21 . 54 - 23.71 

22.60 

2.455 

± 

0.011 

+ 0.046 
-0.082 



2.539 

± 

0.004 

+0.054 

-0.090 



23 . 71 - 26.10 

24.87 

1.883 

± 

0.010 

+ 0.036 
-0.063 



1.972 

± 

0.004 

+0.042 

-0.070 



26 . 10 - 28.73 

27.38 

1.475 

± 

0.008 

+ 0.029 
-0.050 



1.524 

± 

0.003 

+0.033 

-0.054 



28 . 73 - 31.62 

30.13 

1.137 

± 

0.007 

+ 0.022 
-0.038 



1.177 

± 

0.003 

+0.026 

-0.042 



31 . 62 - 35.48 

33.48 

8.629 

± 

0.051 

+ 0.171 
-0.292 

X 

10 "^ 

8.816 

± 

0.019 

+0.195 

-0.318 

X 

10 "^ 

35 . 48 - 39.81 

37.57 

6.254 

± 

0.041 

+ 0.128 
-0.214 

X 

10 ”^ 

6.472 

± 

0.016 

+0.147 

-0.236 

X 

10 "^ 

39 . 81 - 44.67 

42.15 

4.581 

± 

0.033 

+ 0.096 
-0.158 

X 

10 “^ 

4.707 

± 

0.013 

+0.108 

-0.172 

X 

10 "^ 

44 . 67 - 50.12 

47.30 

3.350 

± 

0.027 

+ 0.073 
-0.118 

X 

10 "^ 

3.434 

± 

0.010 

+0.082 

-0.128 

X 

10 "^ 

50 . 12 - 56.23 

53.07 

2.471 

± 

0.022 

+ 0.056 
-0.088 

X 

10 “^ 

2.491 

± 

0.008 

+0.061 

-0.093 

X 

10 "^ 

56 . 23 - 63.10 

59.54 

1.791 

± 

0.017 

+ 0.043 
-0.065 

X 

10 "^ 

1.823 

± 

0.007 

+0.047 

-0.070 

X 

10 "^ 

63 . 10 - 70.79 

66.81 

1.303 

± 

0.014 

+ 0.032 
-0.048 

X 

10 "^ 

1.327 

± 

0.005 

+0.035 

-0.052 

X 

10 "^ 

70 . 79 - 79.43 

74.96 

9.673 

± 

0.114 

+ 0.260 
-0.371 

X 

10 "^ 

9.681 

± 

0.043 

+0.276 

-0.390 

X 

10 "^ 

79 . 43 - 89.13 

84.10 

6.981 

± 

0.092 

+ 0.196 
-0.274 

X 

10 “^ 

7.055 

± 

0.035 

+0.209 

-0.289 

X 

10 "^ 

89 . 13 - 100.0 

94.36 

4.959 

± 

0.073 

+ 0.153 
-0.205 

X 

10 “^ 

5.165 

± 

0.028 

+0.167 

-0.222 

X 

10 "^ 

100 . 0 - 113.6 

106.5 

3.779 

± 

0.059 

+ 0.128 
-0.168 

X 

10 ”^ 

3.736 

± 

0.023 

+0.133 

-0.174 

X 

10 "^ 

113 . 6 - 129.2 

121.1 

2.639 

± 

0.046 

+ 0.101 
-0.125 

X 

10 “^ 

2.605 

± 

0.018 

+0.104 

-0.130 

X 

10 "^ 

129 . 2 - 146.8 

137.6 

1.861 

± 

0.037 

+ 0.080 
-0.096 

X 

10 ”^ 

1.881 

± 

0.014 

+0.084 

-0.101 

X 

10 "^ 

146 . 8 - 166.8 

156.4 

1.322 

± 

0.029 

+ 0.061 
-0.072 

X 

10 “^ 

1.329 

± 

0.011 

+0.064 

-0.075 

X 

10 "^ 
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Table 2 

Helium flux at the top of atmosphere 


BESS-Polar I BESS-Polar II 


Energy Range Ej , Flux ± AFgta ± AFays Flux ± AFsta ± AFgys 


(GeV/nucleon) 

(GeV/nucleon) 

(m^ 

sr 

s GeV/nucleon) 

-1 

(m^ 

sr 

s GeV/nucleon) 

-1 

0.150-0.159 

0.154 

1.346 

± 

0.018 

± 

0.037 

X 

10^ 

3.241 

± 

0.009 

± 

0.082 

X 

10^ 

0.159-0.168 

0.163 

1.406 

± 

0.018 

± 

0.038 

X 

10^ 

3.312 

± 

0.009 

± 

0.084 

X 

10^ 

0.168-0.178 

0.173 

1.470 

± 

0.017 

± 

0.039 

X 

10^ 

3.367 

± 

0.009 

± 

0.086 

X 

10^ 

0.178-0.188 

0.183 

1.540 

± 

0.017 

± 

0.041 

X 

10^ 

3.431 

± 

0.009 

± 

0.087 

X 

10^ 

0.188-0.200 

0.194 

1.601 

± 

0.017 

± 

0.042 

X 

10^ 

3.440 

± 

0.008 

± 

0.088 

X 

10^ 

0.200-0.211 

0.205 

1.583 

± 

0.016 

± 

0.041 

X 

10^ 

3.429 

± 

0.008 

± 

0.088 

X 

10^ 

0.211-0.224 

0.218 

1.662 

± 

0.016 

± 

0.043 

X 

10^ 

3.398 

± 

0.008 

± 

0.088 

X 

10^ 

0.224-0.237 

0.231 

1.654 

± 

0.015 

± 

0.042 

X 

10^ 

3.372 

± 

0.007 

± 

0.087 

X 

10^ 

0.237-0.251 

0.244 

1.674 

± 

0.015 

± 

0.043 

X 

10^ 

3.308 

± 

0.007 

± 

0.086 

X 

10^ 

0.251-0.266 

0.259 

1.706 

± 

0.014 

± 

0.043 

X 

10^ 

3.245 

± 

0.007 

± 

0.085 

X 

10^ 

0.266-0.282 

0.274 

1.676 

± 

0.014 

± 

0.042 

X 

10^ 

3.184 

± 

0.007 

± 

0.083 

X 

10^ 

0.282-0.298 

0.290 

1.675 

± 

0.013 

± 

0.042 

X 

10^ 

3.080 

± 

0.006 

± 

0.081 

X 

10^ 

0.298-0.316 

0.307 

1.646 

± 

0.013 

± 

0.042 

X 

10^ 

3.016 

± 

0.006 

± 

0.079 

X 

10^ 

0.316-0.335 

0.326 

1.630 

± 

0.012 

± 

0.041 

X 

10^ 

2.911 

± 

0.006 

± 

0.076 

X 

10^ 

0.335-0.355 

0.345 

1.605 

± 

0.012 

± 

0.041 

X 

10^ 

2.818 

± 

0.005 

± 

0.074 

X 

10^ 

0.355-0.376 

0.365 

1.604 

± 

0.011 

± 

0.041 

X 

10^ 

2.731 

± 

0.005 

± 

0.072 

X 

10^ 

0.376-0.398 

0.387 

1.573 

± 

0.011 

± 

0.040 

X 

10^ 

2.620 

± 

0.005 

± 

0.069 

X 

10^ 

0.398-0.422 

0.410 

1.531 

± 

0.010 

± 

0.039 

X 

10^ 

2.530 

± 

0.005 

± 

0.067 

X 

10^ 

0.422-0.447 

0.434 

1.498 

± 

0.010 

± 

0.038 

X 

10^ 

2.416 

± 

0.004 

± 

0.065 

X 

10^ 

0.447-0.473 

0.460 

1.467 

± 

0.009 

± 

0.038 

X 

10^ 

2.319 

± 

0.004 

± 

0.062 

X 

10^ 

0.473-0.501 

0.487 

1.425 

± 

0.009 

± 

0.037 

X 

10^ 

2.217 

± 

0.004 

± 

0.060 

X 

10^ 

0.501-0.531 

0.516 

1.382 

± 

0.009 

± 

0.036 

X 

10^ 

2.110 

± 

0.004 

± 

0.058 

X 

10^ 

0.531-0.562 

0.546 

1.338 

± 

0.008 

± 

0.035 

X 

10^ 

2.010 

± 

0.004 

± 

0.056 

X 

10^ 

0.562-0.596 

0.579 

1.272 

± 

0.008 

± 

0.034 

X 

10^ 

1.908 

± 

0.003 

± 

0.053 

X 

10^ 

0.596-0.631 

0.613 

1.230 

± 

0.007 

± 

0.033 

X 

10^ 

1.805 

± 

0.003 

± 

0.051 

X 

10^ 

0.631-0.668 

0.649 

1.195 

± 

0.007 

± 

0.032 

X 

10^ 

1.701 

± 

0.003 

± 

0.049 

X 

10^ 

0.668-0.708 

0.688 

1.133 

± 

0.007 

± 

0.031 

X 

10^ 

1.604 

± 

0.003 

± 

0.047 

X 

10^ 

0.708-0.750 

0.729 

1.079 

± 

0.006 

± 

0.030 

X 

10^ 

1.511 

± 

0.003 

± 

0.044 

X 

10^ 

0.750-0.794 

0.772 

1.026 

± 

0.006 

± 

0.028 

X 

10^ 

1.419 

± 

0.003 

± 

0.042 

X 

10^ 

0.794-0.841 

0.818 

9.866 

± 

0.055 

± 

0.273 

X 

10^ 

1.326 

± 

0.002 

± 

0.039 

X 

10^ 

0.841-0.891 

0.866 

9.291 

± 

0.052 

± 

0.258 

X 

10^ 

1.240 

± 

0.002 

± 

0.037 

X 

10^ 

0.891-0.944 

0.917 

8.909 

± 

0.050 

± 

0.248 

X 

10^ 

1.153 

± 

0.002 

± 

0.034 

X 

10^ 

0.944-1.000 

0.972 

8.162 

± 

0.046 

± 

0.228 

X 

10^ 

1.071 

± 

0.002 

± 

0.032 

X 

10^ 

1.000-1.080 

1.039 

7.661 

± 

0.037 

± 

0.214 

X 

10^ 

9.850 

± 

0.016 

± 

0.292 

X 

10^ 

1.080-1.166 

1.122 

6.971 

± 

0.034 

± 

0.195 

X 

10^ 

8.855 

± 

0.014 

± 

0.262 

X 

10^ 

1.166-1.259 

1.212 

6.444 

± 

0.032 

± 

0.180 

X 

10^ 

7.929 

± 

0.013 

± 

0.234 

X 

10^ 

1.259-1.359 

1.308 

5.773 

± 

0.029 

± 

0.161 

X 

10^ 

7.099 

± 

0.012 

± 

0.210 

X 

10^ 

1.359-1.468 

1.413 

5.209 

± 

0.027 

± 

0.145 

X 

10^ 

6.342 

± 

0.011 

± 

0.187 

X 

10^ 

1.468-1.585 

1.525 

4.724 

± 

0.025 

± 

0.132 

X 

10^ 

5.613 

± 

0.010 

± 

0.166 

X 

10^ 

1.585-1.711 

1.647 

4.188 

± 

0.022 

± 

0.116 

X 

10^ 

4.962 

± 

0.009 

± 

0.146 

X 

10^ 

1.711-1.848 

1.778 

3.749 

± 

0.020 

± 

0.104 

X 

10^ 

4.360 

± 

0.008 

± 

0.128 

X 

10^ 

1.848-1.995 

1.920 

3.331 

± 

0.018 

± 

0.092 

X 

10^ 

3.840 

± 

0.007 

± 

0.113 

X 

10^ 

1.995-2.154 

2.073 

2.944 

± 

0.017 

± 

0.081 

X 

10^ 

3.376 

± 

0.007 

± 

0.099 

X 

10^ 

2.154-2.326 

2.239 

2.616 

± 

0.015 

± 

0.071 

X 

10^ 

2.953 

± 

0.006 

± 

0.086 

X 

10^ 

2.326-2.512 

2.417 

2.274 

± 

0.014 

± 

0.062 

X 

10^ 

2.572 

± 

0.005 

± 

0.075 

X 

10^ 

2.512-2.712 

2.610 

2.004 

± 

0.012 

± 

0.054 

X 

10^ 

2.227 

± 

0.005 

± 

0.065 

X 

10^ 

2.712-2.929 

2.818 

1.759 

± 

0.011 

± 

0.047 

X 

10^ 

1.936 

± 

0.004 

± 

0.056 

X 

10^ 

2.929-3.162 

3.043 

1.522 

± 

0.010 

± 

0.041 

X 

10^ 

1.671 

± 

0.004 

± 

0.048 

X 

10^ 

3.162-3.481 

3.318 

1.297 

± 

0.008 

± 

0.035 

X 

10^ 

1.414 

± 

0.003 

± 

0.041 

X 

10^ 

3.481-3.831 

3.652 

1.077 

± 

0.007 

± 

0.029 

X 

10^ 

1.167 

± 

0.003 

± 

0.034 

X 

10^ 

3.831-4.217 

4.019 

8.975 

± 

0.060 

± 

0.236 



9.528 

± 

0.023 

± 

0.272 



4.217-4.642 

4.424 

7.281 

± 

0.052 

± 

0.191 



7.807 

± 

0.020 

± 

0.223 



4.642-5.109 

4.869 

6.050 

± 

0.045 

± 

0.159 



6.403 

± 

0.018 

± 

0.182 



5.109-5.623 

5.360 

4.949 

± 

0.039 

± 

0.130 



5.190 

± 

0.015 

± 

0.148 



5.623-6.190 

5.899 

3.934 

± 

0.033 

± 

0.103 



4.157 

± 

0.013 

± 

0.118 



6.190-6.813 

6.493 

3.228 

± 

0.029 

± 

0.085 



3.355 

± 

0.011 

± 

0.096 
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Table 2 —Continued 


_ BESS-Polar I _ BESS-Polar II 

Energy Range Ek Flux ± AE^ta ± AF^ys Flux ± AF^ta ± AF^ya 

(GeV/nucleon) (GeV/nucleon) (m^ sr s GeV/nucleon) “ ^ (m^ sr s GeV/nucleon) 


6 . 813 - 7.499 

7.147 

2.574 

± 

0.025 

± 

0.068 


2.700 

± 

0.009 

± 

0.077 



7 . 499 - 8.254 

7.866 

2.044 

± 

0.021 

± 

0.054 


2.159 

± 

0.008 

± 

0.061 



8 . 254 - 9.085 

8.658 

1.675 

± 

0.018 

± 

0.045 


1.727 

± 

0.007 

± 

0.049 



9 . 085 - 10.00 

9.530 

1.288 

± 

0.015 

± 

0.035 


1.355 

± 

0.006 

± 

0.039 



10 . 00 - 11.22 

10.59 

1.020 

± 

0.012 

± 

0.028 


1.061 

± 

0.004 

± 

0.030 



11 . 22 - 12.59 

11.88 

7.776 

± 

0.096 

± 

0.216 X 

10 “^ 

7.912 

± 

0.036 

± 

0.228 

X 

10 

12 . 59 - 14.13 

13.33 

5.714 

± 

0.078 

± 

0.162 X 

10-1 

5.970 

± 

0.030 

± 

0.172 

X 

10 

14 . 13 - 15.85 

14.96 

4.305 

± 

0.064 

± 

0.125 X 

10“1 

4.478 

± 

0.025 

± 

0.130 

X 

10 

15 . 85 - 17.78 

16.78 

3.230 

± 

0.052 

± 

0.096 X 

10“1 

3.330 

± 

0.020 

± 

0.097 

X 

10 

17 . 78 - 19.95 

18.83 

2.415 

± 

0.043 

± 

0.075 X 

10-1 

2.498 

± 

0.016 

± 

0.074 

X 

10 

19 . 95 - 22.39 

21.13 

1.697 

± 

0.034 

± 

0.054 X 

10-1 

1.826 

± 

0.013 

± 

0.055 

X 

10 

22 . 39 - 25.12 

23.71 

1.322 

± 

0.028 

± 

0.044 X 

10-1 

1.377 

± 

0.011 

± 

0.042 

X 

10 

25 . 12 - 28.18 

26.60 

9.943 

± 

0.234 

± 

0.343 X 

10 -^ 

1.004 

± 

0.009 

± 

0.031 

X 

10 

28 . 18 - 31.62 

29.84 

7.177 

± 

0.188 

± 

0.258 X 

10 -^ 

7.565 

± 

0.072 

± 

0.240 

X 

10 

31 . 62 - 38.31 

34.77 

4.778 

± 

0.110 

± 

0.181 X 

10 -^ 

4.901 

± 

0.042 

± 

0.161 

X 

10 

38 . 31 - 46.42 

42.12 

2.837 

± 

0.077 

± 

0.113 X 

10 -^ 

2.950 

± 

0.029 

± 

0.101 

X 

10 

46 . 42 - 56.23 

51.03 

1.750 

± 

0.055 

± 

0.074 X 

10 -^ 

1.774 

± 

0.021 

± 

0.066 

X 

10 

56 . 23 - 68.13 

61.83 

1.121 

± 

0.040 

± 

0.049 X 

10 -^ 

1.109 

± 

0.015 

± 

0.044 

X 

10 

68 . 13 - 82.54 

74.91 

6.627 

± 

0.283 

± 

0.320 X 

10 -^ 

6.546 

± 

0.104 

± 

0.293 

X 

10 
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Table 3 

Proton/Helium ratio at the top of atmosphere 





BESS-Polar I 



BESS-Polar II 



Rigidity Range 

R 

Ratio ± Ai?Osta i AROsys 

Ratio ± AROsta i AROsys 

(GV) 

(GV) 



(one) 




(one) 



1.098-1.132 

1.115 

1.638 

± 

0.022 ± 0.056 

X 10^ 

1.452 

± 

0-004 ± 0.051 

X 

10^ 

1.132-1.168 

1.150 

1.551 

± 

0.020 ± 0.051 

X 10^ 

1-375 

± 

0.004 ± 0.048 

X 

10^ 

1.168-1.205 

1.187 

1.461 

± 

0.018 ± 0.048 

X 10^ 

1.306 

± 

0.003 ± 0.046 

X 

10^ 

1.205-1.244 

1.224 

1.374 

± 

0.015 ± 0.044 

X 10^ 

1-237 

± 

0.003 ± 0.042 

X 

10^ 

1.244-1.283 

1.263 

1.294 

± 

0.014 ± 0.040 

X 10^ 

1.189 

± 

0.003 ± 0.041 

X 

10^ 

1.283-1.325 

1.304 

1.281 

± 

0.013 ± 0.039 

X 10^ 

1.145 

± 

0.003 ± 0.039 

X 

10^ 

1.325-1.367 

1.346 

1.191 

± 

0.012 ± 0.036 

X 10^ 

1-107 

± 

0.003 ± 0.037 

X 

10^ 

1.367-1.412 

1.389 

1.168 

± 

0.011 ± 0.034 

X 10^ 

1.068 

± 

0.002 ± 0.036 

X 

10^ 

1.412-1.458 

1.435 

1.123 

± 

0.010 ± 0.033 

X 10^ 

1.039 

± 

0.002 ± 0.035 

X 

10^ 

1.458-1.506 

1.482 

1.075 

± 

0.009 ± 0.031 

X 10^ 

1.012 

± 

0.002 ± 0.033 

X 

10^ 

1.506-1.555 

1.530 

1.060 

± 

0.009 ± 0.030 

X 10^ 

9-805 

± 

0.021 ± 0.322 



1.555-1.607 

1.581 

1.024 

± 

0.008 ± 0.029 

X 10^ 

9-634 

± 

0.020 ± 0.312 



1.607-1.661 

1.634 

1.008 

± 

0.008 ± 0.028 

X 10^ 

9.326 

± 

0.019 ± 0.298 



1.661-1.717 

1.688 

9.780 

± 

0.075 ± 0.269 


9.133 

± 

0.019 ± 0.292 



1.717-1.774 

1.745 

9.557 

± 

0.072 ± 0.259 


8.915 

± 

0.018 ± 0.282 



1.774-1.835 

1.805 

9.150 

± 

0.066 ± 0.246 


8-670 

± 

0.017 ± 0.272 



1.835-1.898 

1.866 

8.907 

± 

0.063 ± 0.240 


8.510 

± 

0.017 ± 0.267 



1.898-1.964 

1.930 

8.797 

± 

0.061 ± 0.237 


8.295 

± 

0.016 ± 0.260 



1.964-2.032 

1.997 

8.544 

± 

0.058 ± 0.230 


8-147 

± 

0.016 ± 0.256 



2.032-2.103 

2.067 

8.319 

± 

0.056 ± 0.225 


7-961 

± 

0.015 ± 0.250 



2.103-2.178 

2.140 

8.167 

± 

0.053 ± 0.222 


7-793 

± 

0.015 ± 0.247 



2.178-2.255 

2.216 

7.945 

± 

0.051 ± 0.217 


7-657 

± 

0.014 ± 0.245 



2.255-2.336 

2.295 

7.795 

± 

0.049 ± 0.217 


7-505 

± 

0.014 ± 0.244 



2.336-2.421 

2.378 

7.760 

± 

0.049 ± 0.215 


7-376 

± 

0.014 ± 0.238 



2.421-2.510 

2.465 

7.580 

± 

0.047 ± 0.213 


7-267 

± 

0.013 ± 0.239 



2.510-2.601 

2.555 

7.338 

± 

0.045 ± 0.209 


7-171 

± 

0.013 ± 0.239 



2.601-2.699 

2.650 

7.277 

± 

0.044 ± 0.210 


7-050 

± 

0.013 ± 0.238 



2.699-2.800 

2.749 

7.190 

± 

0.043 ± 0.207 


6.933 

± 

0.013 ± 0.234 



2.800-2.907 

2.853 

7.059 

± 

0.043 ± 0.206 


6.836 

± 

0.013 ± 0.233 



2.907-3.018 

2.961 

6.882 

± 

0.041 ± 0.202 


6-751 

± 

0.013 ± 0.232 



3.018-3.134 

3.075 

6.854 

± 

0.041 ± 0.202 


6-664 

± 

0.013 ± 0.230 



3.134-3.256 

3.194 

6.692 

± 

0.040 ± 0.197 


6.602 

± 

0.013 ± 0.228 



3.256-3.385 

3.319 

6.788 

± 

0.041 ± 0.201 


6-547 

± 

0.013 ± 0.226 



3.385-3.566 

3.474 

6.668 

± 

0.035 ± 0.196 


6-473 

± 

0.011 ± 0.223 



3.566-3.759 

3.661 

6.612 

± 

0.035 ±8:}ii 


6-426 

± 

0.011 +8:11? 



3.759-3.966 

3.861 

6.499 

± 

0.034 +g:}i0 


6-397 

± 

0.011 ±8:22? 



3.966-4.186 

4.075 

6.492 

± 

0.035 ±g:}i| 


6-367 

± 

0.011 ±8:224 



4.186-4.423 

4.303 

6.418 

± 

0-035 


6-308 

± 

0.011 ±8:1M 



4.423-4.676 

4.547 

6.316 

± 

0.035 


6.289 

± 

0.012 ±8:lil 



4.676-4.947 

4.810 

6.315 

± 

0.035 


6.242 

± 

0.012 +8:li? 



4.947-5.238 

5.090 

6.247 

± 

0.036 


6-227 

± 

0.012 +8:12f 



5.238-5.550 

5.391 

6.173 

± 

0-036 


6.169 

± 

0.013 +8:lSi 



5.550-5.885 

5.715 

6.152 

± 

0.037 


6.120 

± 

0.013 ±8:182 



5.885-6.244 

6.061 

6.047 

± 

0-037 


6.065 

± 

0.013 ±8:188 



6.244-6.630 

6.434 

6.060 

± 

0-039 ±8iii 


6.033 

± 

0.013 ±8:118 



6.630-7.046 

6.834 

5.968 

± 

0-039 ±8iii 


6.003 

± 

0.014 ±8:lii 



7.046-7.493 

7.265 

5.914 

± 

0-040 ±8iii 


5-943 

± 

0.014 ±8:118 



7.493-7.973 

7.729 

5.866 

± 

0-041 ig:!™ 


5-897 

± 

0.015 ±8:111 



7.973-8.626 

8.292 

5.802 

± 

0-038 


5.832 

± 

0.014 ±8:lii 



8.626-9.342 

8.975 

5.776 

± 

0-040 ±g:Mi 


5-778 

± 

0.014 ±8:lil 



9.342-10.13 

9.725 

5.663 

± 

0-041 ±8:^} 


5-744 

± 

0.015 ±8:11? 



10.13-10.99 

10.55 

5.660 

± 

0-043 


5-688 

± 

0.016 ±8:118 



10.99-11.94 

11.45 

5.540 

± 

0-044 +8:ifi 


5-584 

± 

0.016 ±8:111 



11.94-12.97 

12.44 

5.454 

± 

0-046 


5.515 

± 

0.017 ±8:118 



12.97-14.12 

13.54 

5.505 

± 

0-050 


5.525 

± 

0.018 ±8:11? 



14.12-15.38 

14.73 

5.389 

± 

0-051 


5.441 

± 

0.019 ±8:111 
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Table 3 —Continued 


BESS-Polar I BESS-Polar II 


Rigidity Range R Ratio ± AROsta i AROsys Ratio ± AROsta i AROsys 

(GV) (GV) (one) (one) 


15.38-16.76 

16.05 

5.342 

± 

0.054 

+ 0.142 

5.368 

± 

0.020 

+ 0.170 

16.76-18.27 

17.49 

5.329 

± 

0.058 

+ 0.143 
-0.209 

5.305 

± 

0.021 

+ 0.169 
-0.232 

18.27-19.95 

19.09 

5.208 

± 

0.060 

+ 0.141 

5.225 

± 

0.022 

+ 0.167 

19.95-21.79 

20.84 

5.336 

± 

0.067 

+ 0.146 

5.234 

± 

0.024 

+ 0.168 

21.79-24.23 

22.97 

5.125 

± 

0.063 

+ 0.143 

5.121 

± 

0.023 

+ 0.165 

24.23-26.98 

25.56 

4.987 

± 

0.066 

+ 0.143 
-0.201 

5.122 

± 

0.025 

+ 0.166 
-0.227 

26.98-30.05 

28.46 

5.094 

± 

0.074 

+ 0.150 
-0.209 

5.037 

± 

0.027 

+ 0.165 
-0.224 

30.05-33.51 

31.72 

4.985 

± 

0.079 

+ 0.152 
-0.208 

4.949 

± 

0.029 

+ 0.165 
-0.222 

33.51-37.39 

35.38 

4.962 

± 

0.086 

+ 0.157 
-0.212 

4.909 

± 

0.032 

+ 0.166 
-0.222 

37.39-41.73 

39.48 

4.802 

± 

0.091 

+ 0.159 
-0.210 

4.829 

± 

0.034 

+ 0.166 
-0.221 

41.73-46.60 

44.08 

5.074 

± 

0.108 

+ 0.177 
-0.228 

4.830 

± 

0.038 

+ 0.170 
-0.224 

46.60-52.07 

49.24 

4.808 

± 

0.110 

+ 0.176 
-0.223 

4.721 

± 

0.040 

+ 0.171 
-0.222 

52.07-58.20 

55.03 

4.721 

± 

0.119 

+ 0.181 
-0.225 

4.707 

± 

0.044 

+ 0.174 
-0.224 

58.20-65.08 

61.52 

4.729 

± 

0.133 

+ 0.193 
-0.235 

4.595 

± 

0.047 

+ 0.179 
-0.226 

65.08-78.46 

71.37 

4.772 

± 

0.118 

+ 0.209 
-0.249 

4.680 

± 

0.043 

+ 0.193 
-0.239 

78.46-94.67 

86.08 

4.713 

± 

0.137 

+ 0.222 
-0.260 

4.641 

± 

0.049 

+ 0.206 
-0.248 

94.67-114.3 

103.9 

4.646 

± 

0.157 

+ 0.245 

4.599 

± 

0.058 

+ 0.231 

114.3-138.1 

125.5 

4.361 

± 

0.169 

+ 0.249 

4.369 

± 

0.063 

+ 0.241 

138.1-166.9 

151.6 

4.376 

± 

0.202 

±8:i?e 

4.510 

± 

0.078 

+ 0.296 
-0.325 
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